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A density matrix is defined using coherent states for space-times with apparent 
horizons. Evolving the density matrix in time gives the origin of Hawking radiation. 
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A new theory is expected to take over at Planck distances as 'quantum effects' of gravity 
start dominating. One of the promising approaches to the theory of quantum gravity is 
the theory of Loop Quantum Gravity (LQG), which is by formulation non-perturbative 
and background independent [1-3]. LQG has a well defined kinematical Hilbert space, and 
though the Hamiltonian constraint remains unsolved, the theory allows for a semiclassical 
£> ' sector of the theory. This includes 'coherent states' [4, 5] which are peaked at classical phase 
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space elements. Using these as a starting point, I defined in a series of papers [6-8] coherent 
states for the Schwarzschild space-time, and derived an origin of entropy using quantum 
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mechanical definition of entropy from density matrices. The exact entropy is a function of 
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the graph used to obtain the LQG phase space variables [9]. The zeroeth order term is 
proportional to the area of the horizon signifying a universality of the Bekenstein-Hawking 
entropy. The proportionality constant and the correction terms bring out the details of the 
graph [8]. 

In this paper we take this new way of finding the origin of entropy a step further by 
evolving the spatial slice in time [10], and observing the evolution of the density matrix in the 
process. This state as of now does not satisfy the Hamiltonian constraint, but one is allowed 
to take an arbitrary initial state, or a wavepacket with appropriate properties, representing 
a macroscopic configuration. The evolution discussed in this paper is semiclassical, i.e. no 
attempt is made to use the full Hamiltonian. 

The quasilocal energy (QLE) of an outside observer, defined in [11] is used as the Hamil- 
tonian to evolve the system. As the time clicks in the observers clock, the Hamiltonian 
evolves the coherent state such that the area of the horizon remains the same as predicted 



by classical physics. However, classically forbidden regions become accessible quantum me- 
chanically, and vertices of the graph hidden behind the horizon in one slice emerge outside 
the horizon in the next slice. This gives a net change in area, and the mass deficit is emitted 
from the black hole. This evolution is not unitary, and the quasi-local energy which is used 
to evolve the slice is not mapped to a Hermitian operator. When matter is coupled to the 
gravitational system, a net flux emerges causing a decay of the horizon. 

In section II we introduce the formalism by describing the coherent state, the black hole 
time slice, the apparent horizon equation, and the density matrix. Section III describes the 
time evolution of the system and gives a derivation of the change in entropy. In section 
IV we give a description of a matter current emergent from behind the horizon. Finally in 
the concluding section we include a discussion about the implications of the non-Hermitian 
evolution. 

II. THE COHERENT STATE IN LQG 

For gravity, finding the canonical variables which describe the physical phase space is an 
odd task as there is no unique time. Nevertheless a fiducial time coordinate can be chosen, 
which breaks the manifest diffeomorphism invariance, restored in the Hilbert space of states 
by imposing constraints. 

The constant time slices are described by the intrinsic metric q a b and the extrinsic cur- 
vature K a b (a,b=l,2,3). The theory can be formulated in terms of the square root of the 
metric, the triads e^ defined thus: 

e„e£ = q ab (1) 

where I represents the internal index for the rotation group SO (3) of the tangent space and 
a, b = 1, 2, 3. The internal group is taken to be SU(2), as this is locally isomorphic to SO(3). 
The theory is then defined in terms of the 'spin connection' T^ = e IJK e h jV a &bK and the 
triads. However, a redefinition of the variables in terms of tangent space densitised triads 
E T a and a corresponding gauge connection A ! a where / represents the SU(2) index simplifies 
the quantisation considerably. 

A I a = T I a -PK ah e Ib Ef = I(det e)ej (2) 

(e„ are the usual triads, K a b is the extrinsic curvature, T^ the associated spin connection, 



/3 the one parameter ambiguity which remains named as the Immirzi parameter.) The 
quantisation of the Poisson algebra of these variables is done by smearing the connection 
along one dimensional edges e of length 5 e of a graph T to get holonomies h e (A). The 
triads are smeared in a set of 2-surface decomposition of the three dimensional spatial slice 
to get the corresponding momentum Pf. The algebra is then represented in a kinematic 
'HUbert space', in which the physical constraints have been 'formally' realised [12]. Once 
the phase space variables have been identified, one can write a coherent state for these [4] 
i.e. minimum uncertainty states peaked at classical values of h e , Pf. In analogy with the 
harmonic oscillator coherent states, where the coherent state is a function of the complexified 
phase space element x — ip, the SU(2) coherent states are peaked at the complexified phase 
space element g e = e tT p ^ 2 h e . These g e are thus elements in the complexification of SU(2) 
as e tT p I 2 (T 1 being the generator matrices of SU(2)) is a Hermitian matrix and h e is the 
unitary SU(2) matrix. Whether these are physical coherent states, or have appropriate 
behavior under the action of the constraints has to be examined carefully [13]. The coherent 
state in the momentum representation for one edge is defined to be 

W{9e) >= E e-M+WicMlJmn > (3) 

jmn 

In the above g e is a complexified classical phase space element e lT p <= c / 2 hf, (the P/ cl and 
the hf represent classical momenta and holonomy obtained by embedding the edge in the 
classical metric). The \jmn > are the usual basis spin network states given by 7ij(h) mn , 
which is the jth representation of the SU(2) element h e . Similarly, (2j + 1) x (2j + 1) 
dimensional representations of the 2x2 matrix g e are denoted as ftj{g e ) mn - The j is the 
quantum number of the SU(2) Casimir operator in that representation, and m,n represent 
azimuthal quantum numbers which run from —j..j. The coherent state is precisely peaked 
with maximum probability at the hf for the variable h e as well as the classical momentum 
pic\ f or ^g variable Pf. The fluctuations about the classical value are controlled by the 
parameter t (the semiclassicality parameter). This parameter is given by I 2 /a where l p is 
Planck's constant and a a dimensional constant which characterises the system. The coherent 
state for an entire slice can be obtained by taking the tensor product of the coherent state 
for each edge which form a graph T, 

*r = IK- (4) 



In [7] the g e was evaluated for the Schwarzschild black hole by embedding a graph on a 
spatial slice with zero intrinsic curvature. The particular graph which was used had the 
edges along the coordinate lines of a sphere. This simplistic graph, was very useful in 
obtaining the description of the space-time in terms of discretised holonomy and momenta. 
A particularly interesting consequence of this was that the phase space variables were finite 
and well defined even at the singularity. 

Given that the area of a surface in gravity is measured as the integral of the square root 



of the metric over the surface, the area operator can be written simply as A = y P^PJ. The 
expectation value of the area operator in the coherent state emerges as [9] 

<^|i|^>=( J + I)t (5) 

Thus we are considering a semiclassical state, which is a state such that expectation values 
of operators are closest to their classical values. The information of the classical phase space 
variables are encoded in the complexified SU(2) elements labeled as g e . The fluctuations over 
the classical values are controlled by the semiclassical parameter t. 

The density matrix which describes the entire black hole slice is obtained as 

p Total = |^ r >< ^ r | ( 6 ) 

where |\l/r > is the coherent state wavefunction for the entire slice, a tensor product of 
coherent state for each edge. 

A. Apparent Horizons 

We concentrate on the coherent state near the apparent horizon contained in the spatial 
slice. We find that motivated from the apparent horizon equation the graph across the 
horizon can be taken to be populated by radial edges, linking vertices outside and inside 
the horizon. One then traces over the coherent state within the horizon. Initially we take 
a particular time slicing of the black hole, which has the spatial slices with zero intrinsic 
curvature [7]. One such metric which has the time slices as flat is the Lemaitre metric 
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r 2/3 (d9 2 + sin 2 9d<j) 2 ). (7) 



The r g = 2GM, (in units of c=l)and in the r = constant slices one can define the induced 
metric in terms of a V coordinate defined as dr = dRj [3/2r g (R — r c )] (r = r c ) on the 
slice. One gets the metric of the three slice to be 

dsl = dr 2 + r 2 (de 2 + sin 2 9d(f) 2 ). (8) 

The entire curvature of the space-time metric is contained in the extrinsic curvature or 
Kfj,v = \d T g^u tensor of the r = constant slices. Now if there exists an apparent horizon 
somewhere in the above spatial slice, then that is located as a solution to the equation 

V a S a + K ab S a S b -K = (9) 

where S a , ((a, b = 1, 2, 3) denote the spatial indices) is the normal to the horizon, K ab the 
extrinsic curvature in the induced coordinates of the slice, and K the trace of the extrinsic 
curvature. If the horizon is chosen to be the 2-sphere, then in the coordinates of (8), 
S a = (1, 0, 0), the apparent horizon equation as a function of the metric reduces to: 

K rr (l - q rr ) - K m H ~ Tj r - K ee q ee -T e gr = (10) 

Note that the first term of the equation disappears trivially as 1 = q rr for any point in the 
spatial slice. Even at the operator level the q rr can be set to the identity operator in the first 
approximation, as q rr = P Sr P er /V 2 (V being the volume operator) upto normalisations, and 
in the spherically symmetric metric V = P er 8 er (upto discretisation constants). Thus the 
operators in the numerator and denominator cancel and the normalisation conspire, leaving 
q rr = X. To understand the rest of the equation in terms of the holonomy and momentum 
variables of LQG, which are classically measured in the same metric as (8), we use the 
following regularisation 

Ke{4,),e{4>) = ee(<t>) K e(<t>y Qe(<t>)0(<P) = e 9(<p) e e(<t>) ( ll ) 

e I $tW = NTr[T I h- 1 w {h e , w ,V}] (12) 

(N is a constant, a function of the edge lengths and the area bits of the discretisation) 
and V is the volume operator. 

Kfa) = ^^T'^K] (13) 



Here (5 has been used as a parameter to identify the K T a operator, and this is mainly a trick. 
In the continuum limit 

h e (Al) = Limit^oe-f Aadxa = (I + AlT'SeJ (14) 

As the gauge connection is a function of the Immirzi parameter due to (2), the expectation 
value of this operator in a coherent state will be a function of the Immirzi parameter. By 
taking the derivative wrt to the Immirzi parameter we are giving the same status to the 
parameter as is given to 'dimension' in a dimensional regularisation of Feynman diagrams. 
We let the parameter vary by an infinitesimal amount from its value in the particular quan- 
tisation sector, take the derivative, and put its original value in the final answer for the K^ 
operator. The formula (13) is facilitated by the fact that the dependence of A T a on the /3 
is linear. One way to check whether this gives the proper answer is to take a solved quan- 
tum mechanical system and use a similar method there. The most useful example is the 
Harmonic Oscillator Hamiltonian, which can be written as 

v 2 1 
H = — + -m u 2 x 2 (15) 

2m 2 K J 

The ground state is a coherent state, so we take that as an example. We define the 

operator 

Thus 

2 2 dH 2d fhu\ h 

< x 2 >= - < — >= -—; — = (17) 

m ouj m ou z \ 2 J 2mu 

The regularisation (13) is thus an allowed approximation. 

The terms involving the Christoffel connections like T e dr include derivatives in the regu- 
larised version, the derivatives appear as difference of triads across two vertices. Thus 

C = e?e?-L (el{ Vl ) - efe)) e J e ( Vl ) (18) 

As a result of this if we impose restrictions on the Christoffel connections and one of the 
vertices v\ is within the horizon, whereas i>2 is outside the horizon, there will be correlations 
across the horizon. 

If one evaluates the expectation value of the apparent horizon equation using the regu- 
larised variables in the coherent states, then one would obtain 



4 < i>\Pe e 



Tr (T J KX /2 K) Vi ~ Tr (^V^X) J Tr (^KX^o)^ IV> > 



~N' < ^|Tr (h^T 1 ^ h e )j P/ e Vl \^ >= (19) 

(N' is a constant) 

B. Density Matrix 

The density matrix is obtained as 

p Total = |^ r >< ^| ^ 20 ) 

where |\l/r > is the coherent state wavefunction for the entire slice, a tensor product of 
coherent state for each edge. 

But given this, we concentrate in a 'local' region to see the behavior of the horizon 

Total outside local inside (01\ 

P — P P P \ zl ) 

where p local covers a band of vertices surrounding the horizon one set on a sphere at radius 
r g ~ <^e r /2 and one set on a sphere at radius r g + S £r /2 within the horizon, as described in 
[9], and in the figure enclosed. This local density matrix and the correlations due to the 
apparent horizon equation (19) was used to derive entropy [6]. This entropy counts the 
number of ways to induce the horizon area using the spin networks, though the constraints 
have not been appropriately imposed as was obtained using a Chern-Simons theory in [14]. 
However, the entropy calculation using the coherent states provides a tracing mechanism, 
and a method to obtain correlations across the horizon which are gravitational in origin. We 
will henceforth deal with p local ; but we will drop the local label for brevity. 



Horizon 



V 



1 radial 




e' 

9 


edge 


e » % 2 



III. TIME EVOLUTION 

In physical systems, the Hamiltonian generates time evolution, but in General Theory 
of Relativity, the Hamiltonian is a constraint and generates diffeomorphisms in the time 
direction. So the question is, what is physical time, and if that exists, what would be 
the operator evolving the system in that direction? In case of space-times with time like 
Killing vectors, notion of time can be identified with the Killing direction, and a notion of 
'quasilocal energy' (QLE) defined using the same. The QLE then generates translations 
in the Killing time. In case of the Schwarzschild space-time, the QLE has been defined in 
[11]. We build the Hamiltonian which evolves the horizon from one time slice to the next by 
appropriately regularising the QLE. Note the 'Killing time' and QLE are classical concepts, 
and thus regularising QLE gives us a 'semiclassical' Hamiltonian. 

A. Change in Entropy 



Before we get into the analysis of what QLE evolution means, we take a simple system 
made up of two subsystems, and examine the consequences of a Hamiltonian evolution. Let 
the density matrix be defined for a system whose states are given in the tensor product 



Hilbert space H i ® H 2 and given by 



\i)>=Y J d lj \i> \j > 



(22) 



i.i 



where \i > is the basis in Hi and \j > is the basis in H2 and dij are the non-factorisable 
coefficients of the wavefunction in this basis. Let us label the wavefunction at time t — to 
be given by the coefficients d?-. The density matrix is 



P °= £ 4*4,1* >|j><j'|<*'| 

iji'j' 

The reduced density matrix if one traces over Hi is: 



(23) 



V = EEW><i' 



(24) 
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We now evolve the system using a Hamiltonian which has the matrix elements H^^i > 
\j >< j'\ < i'\, we assume that the Hamiltonian does not factorise, that is there exists 
interaction terms between the two Hilbert spaces. The evolution equation is: 



ih?g = [H, P ] 



(25) 



which in this particular basis gives the density matrix elements at a infinitesimally nearby 
slice to be 



j5t*iSt jO* ^0 r 

a i'j' a ij ~ a i'j' a ij + 0T 



Z^V^ijkldkfii'j' d i: jd kl Hj 



kli'j' 



L kl 



(26) 



Thus we evolve the 'unreduced' density matrix and then trace over the Hi in the evolved 
slice. The reduced density matrix in the evolved slice is: 



V^ jSt*jSt _ V^ rf)* jO 



JOT* JOT 

^ a i'j a ij 



M*M 
/ j u i'j u ij 



i* 



7 , [Hijkid'kfli'j d, i jd kl Hi,i i ij 
mj 



(27) 



This gives: 



P 



p° - ^StA 
n 



(28) 



where A represents the commutator. Clearly the entropy in the evolved slice evaluated 
as Sft U — —Tr(plnp) can be found as 



S& = S° BU + -5T[TrA\np° + Trp p° ~ X A) 



(29) 



10 
Given the definition of An, one gets 

An = Y^ [Pijki H kiij - H ijk ip° klij \ (30) 

jki 

In case both the Hamiltonian and the density operator are Hermitian, one obtains 

J2 Ajj = 2l ImTr(p°H) (31) 

3 

This is clearly calculable, and gives the change in entropy ASbh- The hip term yields 
corrections, and we ignore it in the first approximation. 

B. The Hamiltonian 

To trace the origin of Horizon fluctuations, we must take an observer who is stationed 
outside the horizon, or in other words is not a freely falling observer. The quasilocal energy 
is defined using a 'surface' integral of the extrinsic curvature with which the surface is 
embedded in three space. In our case, we take the bounding surface to be the horizon and 
the quasilocal energy is given by the surface term[ll, 15]. 

H = - f d 2 x^k (32) 

K J 

where k is the extrinsic curvature with which the 2-surface, which in this case is the horizon 
S 2 is embedded in the spatial 3-slices, and a is the determinant of the two metric a^ v defined 
on the 2-surface. This 'quasilocal energy' is measured with reference to a background metric. 
Thus H = H — H . We concentrate on the physics observed in a observer stationed at a 
r = constant sphere. 

The metric in static r = const observer's frame is 

ds 2 = -fdt 2 + r 2 (d6 2 + sin 2 Qd<f). (33) 



The / = J 1 — r g /r where r g is the Schwarzschild radius. If we take n^ to be the space-like 
vector, normal to the 2-surface, then the extrinsic curvature is given by: 

K» = <V Q n„ (34) 

and the trace is obviously 

k = V a n a (35) 
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In the special slicing of the of the stationary observer the normal to the horizon 2-surface 
is given by (0, f(r), 0, 0). However, we built the coherent state on the Lemaitre slice. The 
Lemaitre and the Schwarzschild observer's coordinates are related by the following coordi- 
nate transformations, 

[Ldr = (dR - dr) dt = —^— (dr - f'dR) f = ^ (36) 

The r= const cylinder of the Schwarzschild coordinate corresponds to dR = dr of the 
Lemaitre coordinates, and for these dt = dr. Thus unit translation in the t coordinate 
coincides with unit translation in the r coordinate. Further, the intersection of the r = 
constant cylinder with a t = constant surface coincides with the intersection of r=constant 
and the r=constant surface. Thus in the initial slice, the QLE Hamiltonian can be written 
as 

H=±f de d(t)V m^[-9 9ed -^ - 9**^]f(r) - H (37) 

The reference frames' quasilocal energy is a number, it just defines the zero point Hamilto- 
nian. Thus, we replace the classical expressions by operators evaluated at the r = constant 



slice. In the first approximation we simply take the f(r) as classical Jl — r g /r = JS er /2r g = 
e, as this arises due to the coordinate transformation and the norm of the vector n r in the 
previous frame. In the re-writing of (37) in regularised LQG variables the Hamiltonian 
appears rather complicated. 

One can rewrite these in a much simpler form, using the apparent horizon equation. 
Since the Hamiltonian is an integral over the horizon, the variables will satisfy the apparent 
horizon Equation (10) upto quantum fluctuations. Thus the Hamiltonian operator is then 
re-written as 

^horizon = - / d6d<P^9TA K y e + K l^\ ( 38 ) 

where we have used the classical apparent horizon equation (10) (with q rr = 1). 

^Horizon = ^^~ £ Tr^T 7 /?— h ee \P egl + h.c. + (6 -> 0) (39) 

ZKO eg S ee Vl Op 

where C consists of some dimensionless constants s e0 is the 2-dimensional area bit over which 
Ej is smeared, a is a dimensionfull constant which appears to get the P eg i dimension less. 
5 eg is the length for the angular edge eg over which the gauge connection is integrated to 
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obtain the holonomy. The sum over v\ is the set of vertices immediately outside the horizon. 
The (39) can then be lifted to an operator. 

This regularised expression for QLE is for the horizon 2-surface only and would not apply 
for any other spherical surface in the Schwarzs child space-time. 

C. U(l) Case 

Let us take the U(l) case to make the calculations easier and observe the action of 
the QLE Hamiltonian on the evolution of the coherent state. The spin network states are 
replaced by \n >= e m< > , < ( < 2n, n is an integer and the coherent states are: 

^\g e ) = ]T e -(*™ 2 )/2 e ™(xe-ip e ) e -mC ( 40 ) 

n 

g n e = e m<yXe ~ vpe ' > is the complexified phase space element in the 'nth' representation. 
The QLE operator also takes the simplified form 

1 f) 1 f)h~^ 

HS£L = -^C'iK l ^—h e p e + -CTiM-Jj-he (41) 

The prefactors have been clubbed into C . 

In the calculation of the matrix elements, we drop the label of the edges e for the Hamil- 
tonian. 

< m\K^on\n >= J e-^HjgLe^dC (42) 

This calculation can be done by putting an assumption that the C = Ci + $(2- I n this 
Ci,C2 are completely independent of (5. It is an allowed assumption, and identifies the /3 
dependence of the operator matrix elements, which are otherwise 'hidden'. The calculation 
however introduces an arbitrariness in the formula, which can be fixed by requiring that the 
expectation value of the Hamiltonian agrees with the classical QLE [10]. However, in this 
paper we use the 'annihilation' operators defined in [16]. 

This is done by observing that the U(l) coherent states are eigenstates of an annihilation 
operator defined thus: 

g e = e 4/ V e h e g e \i> >= ge\*P > (43) 

The holonomy operator can thus be written as 

k = e-^e'^ g e (44) 
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And the derivative wrt Immirzi parameter of the holonomy which appears in the definition 
of the Hamiltonian replaced by 



fl- 



ap 



e -t/2 



-t/2 






(45) 
(46) 



The dependence of the operator p on the Immirzi parameter is known (2), and thus we 
could evaluate the derivative (fidp p e {/3) = f3dp(p e (l) / (3) = —p e {P)) 
The term 

Tr(p°i/ZL) (47) 

is then computable. Let us take the first term of (41) and find (47). As p° = \ip >< ip\, 
(47) gives simply (we drop the 'e' label for brevity) 



r u(i) 



1 m.t-io d 



< mZLM> > = < i>\ - -Cth-'p—h p\^> + < i>\h.c.\i> > 



d/3 



{At 



p e~*g + e-^0 



2 y 



< ip\e p p e p p \ip > g+ < ip\e p fijr^V \i> > 



p\ip > + < ip\h.c.\ip > 

+ < ip\h.c.\ip > 



-2p o d 9 



We then concentrate on the 2nd term of the above 



-2p p dg 



< i>\e~*P-±v\i> > 



an' 
Mg'W-^ < i>\e~ 2p W >< v'IpIV > 



(49) 

(50) 

(51) 
(52) 



Where we have used the fact that coherent states resolve unity. It can be shown that the 
expectation value of the operators in the t — > collapses the integral to g' — g point [16]. 
Thus one obtains from the above 



.Orrtf(l) 



11 \P -"Horizon 



1 



iC'e' 1 ' 2 



p + g*e- 2p (5 



dg_ 

d(3 



p + h.c. 



'fa 



(53) 
(54) 
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which is real, and thus 

ASbh = (55) 

(this is actually the classical QLE as it should be from Tr(p°ff Horizon))- 
This is obvious, as the way the Hamiltonian is defined, this is simply a function of the 
Hilbert space outside the horizon, and the matrix elements of this will not yield anything 
new. We approximated the horizon sphere by summing over v\ vertices immediately outside 
the horizon. We could do the same by summing over t> 2 vertices immediately within the 
horizon. For the Lemaitre slice, the metric is smooth at the horizon, and one can take the 
'quantum operators' evaluated at the vertex V2- In this case however, as the region is within 
the classical horizon, the norm of the Killing vector is negative, and n r has components 
which are imaginary. The e — > ±ie. Thus -^Horizon = ^E^ H V1 +J2v 2 H V2 \. In the evaluation 
of the QLE, the energy would emerge correct in the 5 er — >■ limit as e — > The regularised 
Hamiltonian is not Hermitian, and the evolution equation is 

ih^- = Hp- P H ] (56) 

And thus the operator which appears in the change of entropy equation is 

AS B „ = ^Tr[Hp° - p°H<\ (57) 

As — 4 c '4 p (58) 

The 'rate of change' of entropy is thus 

AS BH = ^ P (59) 

(we extracted the k from C to get ti and rewrote the rest of the constants as C). 

Thus there is a net change in entropy, but, to see if this is Hawking radiation, we have 
to couple matter to the system. 

D. SU(2) Case 

The SU(2) case is easily reduced to the U(l) case in the actual calculation due to the 
gauge fixing. This is achieved by making the following observations: To retain the metric 



15 
as in the same form as the classical metric, we impose the conditions at the operator level 

Pe a -Pe b = (60) 

such that the corresponding metric has only the diagonal terms as non-zero. With these 
additional 'constraints' on the operators, we can put the P/ such that each has only one 
component surviving, let's say P/ = S^P eg . This also makes the holonomy restricted to the 
U(l) case, as the gauge connection A\ gets restricted to the 1 = 3 and other directions can 
be put to zero. Thus we can take the holonomy to be diagonal 



he = I (61) 




The operator is then obtained as 

H = TrK^T^^P/ = Tr[/,-T 3 /3^]P e 3 = /3^P e 3 (62) 

This is same as the U(l) Hamiltonian (upto normalisations). The spin network states also 
project on to U(l) subgroup, thus giving us the same techniques to use in the calculation 
of the U(l) states as for this one. To observe this, the non-zero elements for the holonomy 

matrix 

I a b\ 
h=\ (63) 

\—b a J 

in the jth representation is given by: 

n(h) - V V^ - "W + m)\{j - n)\(j + n)! m 

Clearly in the particular case we are considering, the 6 = 0, and m,n=-j and j. Thus the 
two non-zero elements are 

7r,(%, = e 2 ^ 7r,-(/0_, w = e" 2 ^ (65) 

The sum over j in the Tr(p if H orizon) with the coherent state defined in (3) thus reduces 
to the U(l) case in the computation of the change in entropy. Thus the rate of change in 
entropy of a classically spherically symmetric black hole is given by 



p v 



d Xe p , 9 X e, 

d/3 ee d/3 e * 



(66) 
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where the classical holonomies h ee(4>) = e e w . If we plug in the actual values, we get 
this to be 

AS BH = ±^ 1 E dA ^ r 9 (67) 

p V 

where dA v the area element at vertex v on the sphere. This change in entropy is totally 
gravitational in origin, and seems to signify the emergence of 'geometry' from within the 
horizon. 

In fact, if we some over the area, we get the ASbh = ± 87r ,2 T r„ (if we set C = 1//3), which 
would be the change in entropy when the radius of the horizon changes by 5r g = eSr\ 

IV. OUTGOING FLUX OF RADIATION 

In the previous section we found that as the system evolves in time, the horizon fluctuates 
and the area decreases. But is this Hawking radiation? Adding matter to a 'coherent state' 
description of semiclassical gravity has been discussed [17]. Thus, given a massless scalar 
field Lagrangian coupled to gravity, whose Hamiltonian is given by 

H sc = fd 3 x[^ + (V<P) 2 }, (68) 

the 'gravity' in the Hamiltonian can be regularised in terms of the h e , P/ operators in the 
coherent state formalism. The integral over the three volume gets converted to a sum over 
the vertices dotting the region. Thus 

H: c = J2H v (h e ,P*,V) (69) 

V 

This Hamiltonian is an operator, and one evaluates an expectation value of the Hamil- 
tonian in the reduced density matrix of the initial slice, to find the classical behavior of the 
scalar field as observed by an observer outside the horizon. Thus 

Tr(p^ c ) (70) 

This Hamiltonian and the density matrix are then both evolved according to the time-like 
observers frame. One gets 

% h^ = [H,H sc ] (71) 

This gives 

Tr ( P t+5t H^) - Tr tflTJ = -(5t) 2 Tt{[H, p^][H, H%} (72) 
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It is very clear thus that the order St terms are zero for this. However, allowing for the 
non-unitary evolution using the non- Hermit ian Hamiltonian, the St terms survive. In fact 
the terms are 

- L -^Ti[(Hp - P H^)H SC ] - ^Tr[p(HH sc - H SC H^)} (73) 

The first term is remarkable, it shows that the term giving rise to entropy change teams up 
with the expectation value of the scalar Hamiltonian. The second term yields corrections, 
so we ignore it in the first approximation. The exact details of the computation have to 
be obtained using the coherent state of the matter and gravity coupled system [17]. If one 
simple takes the matter + gravity system in a tensor product form, and one has matter 
quanta of energy u sitting at one vertex, then the first term would give new matter in the 
evolved slice as ASbh w. The 'rate' of particle creation thus has the form —2coo/Th where 
Tjj is the Hawking temperature for the signs +(— )ie and negative (positive) u. 

Thus from the above it seems 
(i)One has found emission of matter quanta from a black hole but from a 'semiclassical' 
description rooted in a theory of quantum gravity, beyond quantum fields in curved space- 
time. 

(ii)The results indicate a non-unitary evolution which allows space to emerge from within 
the horizon. 

(iii) The emission is perceived by a static or an accelerating observer as anticipated, and the 
non-unitary flow might be due to the semiclassical approximations. A quantum evolution 
using the quantum Hamiltonian might still be unitary. 

The above derivation seems to be a 'quantum gravity' description of the tunneling mecha- 
nism for describing Hawking radiation [18]. However, the results are preliminary and further 
investigation has to be done. 

V. CONCLUSION 

In this paper we showed a method to obtain the origin of Hawking radiation using a co- 
herent state description of a black hole space-time. We took a quantum wavefunction defined 
on an initial slice, peaked with maximum probability at classical phase space-variables. We 
then evolved the slice using a Hamiltonian, which is the 'quasilocal energy' at the horizon. 
This QLE evolved the system in the time and the entropy was shown to change, indicating 
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a change in black hole mass and hence an emergence of interesting non-unitary dynamics. 
One of course has to investigate further to see what is the endpoint of this time evolution. 
The time flow indicates one might have to formulate quantum theory of gravity rooted in 
irreversible physics. The presence of additional degrees of freedom in the form of 'graphs' 
also indicates that the classical phase space might not be described by deterministic physics, 
but by distributions, a manifestation of microscopic irreversible physics in complex systems. 
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